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Abstract 

The  morphology  diagram  of  possible  structures  in  two-dimensional  diffusional  growth  is  given  in  the  parameter  space  of 
undercooling  A versus  anisotropy  of  surface  tension  e.  The  building  block  of  the  dendritic  structure  is  a dendrite  with  a 
parabolic  tip,  and  the  basic  element  of  the  seaweed  structure  is  a doublon.  The  transition  between  these  structures  shows  a 
jump  in  the  growth  velocity.  We  show  the  analogy  of  diffusional  growth  with  dewetting  patterns  of  a fluid  film  on  a 
substrate.  We  also  describe  the  structures  and  velocities  of  fractal  dendrites  and  doublons  destroyed  by  noise.  The  extension 
of  these  results  to  three-dimensional  growth  is  briefly  discussed.  © 2000  Elsevier  Science  B.V.  All  rights  reserved. 


1.  Introduction 

During  the  last  few  years,  our  understanding  of 
pattern  formation  in  various  non-linear  dissipative 
systems  has  made  remarkable  progress.  Building  on 
these  foundations,  it  has  now  become  possible  to 
develop  a description  of  a large  class  of  patterns  that 
are  found  in  diffusional  growth.  This  leads  to  the 
construction  of  a morphology  diagram  and  to  predic- 
tions concerning  the  transitions  between  the  different 
structures  [1,2], 

To  introduce  the  topic,  let  us  take  a look  at  a few 
interesting  patterns  appearing  in  nature  and  experi- 
ments. Fig.  1 displays  one  of  the  most  popular 
examples  which  has  become  a paradigm  of  the  field: 
the  snowflake.  It  also  happens  to  be  the  first  case  of 
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Fig.  1.  Drawing  of  a snowflake,  adapted  from  a photograph  by 
Kobayashi  and  Kuroda  [3], 
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spontaneous  pattern  formation  in  nature  that  has  been 
treated  scientifically  (rather  than  as  a theological 
topic)  in  a publication.  In  his  article  on  ‘The  Six- 
Cornered  Snowflake’  [4],  Kepler  speculated  that  the 
six-fold  symmetry  of  snowflakes  has  to  do  with 
something  we  would  call  the  crystal  structure  today. 
He  also  arrived  at  the  conclusion  that  a detailed 
understanding  of  snowflakes  was  beyond  his  reach. 
Just  how  judicious  this  modest  statement  was,  may 
be  gathered  from  the  fact  that  even  10  years  ago, 
more  than  350  years  after  Kepler,  the  basic  mecha- 
nism by  which  the  characteristic  length  scales  of 
snowflake  patterns  arise  was  not  known. 

Snowflakes  are  ice  crystals  growing  from  a critical 
nucleus  in  an  environment  of  supercooled  vapor.  The 
condensation  process  clearly  is  a first-order  phase 
transition.  Contrary  to  second-order  transitions,  first- 
order  ones  give  rise  to  a sharp  interface  between  the 
two  phases.  It  is  the  dynamics  of  this  interface  that 
are  of  interest  here. 

Snowflakes  are  flat  objects  of  approximate  but  not 
exact  hexagonal  symmetry.  Their  six  main  arms  have 
a characteristic  structure:  they  emit  sidebranches  at 
crystallographic  angles.  When  these  structures  were 
given  names,  classical  education  was  still  prevailing, 
so  they  were  called  dendrites,  after  the  greek  'to 
8ev8pov' , meaning  ‘the  tree’. 

No  two  snowflakes  are  alike.  Therefore,  their 
shape  must  be  extremely  sensitive  to  the  environ- 
mental changes  brought  about  by  their  turbulent 
motion  within  the  cloud.  On  the  other  hand,  the  six 
arms  are  very  similar  to  each  other.  Hence,  on  the 
length  scale  of  the  size  of  a snowflake  (1  cm),  its 
environment  must  be  spatially  homogeneous. 

The  question  then  arises  why  such  a complex 
morphology  evolves  in  a uniform  environment  — 
why  not  simply  spheres?  Or,  if  the  crystalline 
anisotropy  plays  a role,  why  not  just  hexagonal 
plates  (according  to  the  crystal  system  of  ice)? 

Another  question  that  we  may  naturally  ask  is: 
how  does  the  snowflake  select  its  typical  length 
scales,  e.g.  the  tip  radius  of  its  main  arms,  the  width 
and  spacing  of  the  sidebranches? 

The  growth  dynamics  of  true  snowflakes  are 
complicated  by  several  circumstances  — there  is  a 
large  density  difference  between  the  solid  and  the 
vapor;  besides  the  vapor  there  is  a second  agent  for 
heat  transport,  namely  air;  furthermore,  the  temporal 
fluctuations  of  the  environment  of  a flake  are  strong. 


Clearly,  it  will  be  much  easier  to  understand  the 
essential  mechanisms  by  analysing  experiments 
which  are  simple  to  prepare  and  still  produce  den- 
drites. Some  most  careful  and  beautiful  experiments 
of  this  kind  have  been  done  in  the  group  of  Glicks- 
man  [5-8].  They  consist  of  chemically  pure  suc- 
cinonitrile  [C2H4(CN)2],  which  is  an  important 
model  substance  because,  being  transparent,  it  can  be 
easily  observed  during  solidification,  while  its  solidi- 
fication properties  are  similar  to  those  of  (technologi- 
cally important)  metals.  The  solid  phase  grows  into 
an  undercooled  melt.  There  are  no  facets  on  the 
crystal,  its  surface  is  rough.  This  means  that  the 
attachment  kinetics  of  molecules  at  the  interface  are 
very  fast  — the  free  energy  for  the  formation  of  a 
step  on  the  interface  is  zero,  there  is  no  nucleation 
banner.  Thus  kinetics  are  not  a limiting  factor  in  the 
growth  process.  Growth  is  controlled  by  how  fast 
latent  heat  produced  in  the  solidification  process  can 
be  transported  away  via  diffusion.  (In  the  case  of  the 
snowflake,  this  is  not  true  for  the  third  dimension; 
growth  in  that  direction  is  faceted,  keeping  the  flake 
essentially  two-dimensional.) 

The  main  result  of  these  experiments  to  be  kept  in 
mind  is  that,  given  the  undercooling,  the  dendrite  has 
both  a uniquely  selected  tip  radius  and  constant 
velocity. 

Completely  different  structures  are  commonly 
formed  by  minerals  crystallizing  from  viscous  mag- 
mas or  by  certain  polymers  solidifying  from  the  melt 

[9] ,  so-called  spherulites.  Contrary  to  dendrites, 
spherulites  are  polycrystalline.  They  have  approxi- 
mately radial  symmetry  and  their  substructures  are 
not  oriented  along  crystallographic  directions. 

Their  radius  grows  proportional  with  time,  which 
means  they  have  a stationary  growth  phase.  Again, 
one  would  like  to  be  able  to  make  predictions  about 
selected  length  scales. 

It  was  proposed  by  Goldenfeld  [9]  that  spherulites 
are  three-dimensional  examples  for  a growth  struc- 
ture which  was  then  called  ‘dense  branching  mor- 
phology’ and  will  be  described  shortly.  Experimen- 
tally, densely  branching  patterns  were  first  obtained 
in  a non-crystalline  system;  this  was  the  viscous 
fingering  experiment  in  a circular  Hele-Shaw  cell 

[10] .  An  inviscid  fluid  is  injected  into  a viscous  one, 
displacing  it  and  forming  branchy  structures  in  the 
process.  It  is  still  controversial  [11,12]  whether  in 
such  a system,  governed  by  the  Laplace  instead  of 
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the  diffusion  equation,  the  dense  branching  morphol- 
ogy can  persist  at  large  length  scales. 

There  are  well-known  growth  models,  most  no- 
tably diffusion-limited  aggregation  (DLA),  that  lead 
to  the  formation  of  fractal  patterns  [13], 

Since  DLA  is  a Laplacian  system,  it  is  natural  to 
ask  whether  fractal  structures  are  to  be  expected  in 
diffusional  growth,  too. 

An  obvious  classification  of  growth  structures 
would  then  be  to  distinguish  between  compact  and 
fractal  patterns  [1,2].  By  compact  growth  we  mean 
growth  at  a constant  (average)  density,  irrespective  of 
the  value  of  this  density.  As  we  shall  see  later,  fractal 
patterns  in  crystal  growth  are  fractal  only  up  to  a 
certain  size  and  compact  beyond. 

What  we  are  aiming  at  is  to  formulate  a theory 
that  relates  the  diverse  discussed  patterns  in  some 
kind  of  kinetic  phase  diagram.  The  ‘kinetic’  used 
here  reminds  us  that  we  are  not  dealing  with  phase 
transitions  (which  would  imply  thermodynamic 
equilibrium),  but  dynamic  states  manifesting  them- 
selves in  typical  growth  patterns.  In  the  literature,  the 
term  morphology  diagram  is  widely  used  [14], 

On  the  one  hand,  such  a diagram  indicates  which 
patterns  are  selected  for  which  system  parameters.  Its 
existence  signifies,  in  particular,  that  there  is  a 
certain  degree  of  independence  of  the  observed 
morphologies  of  initial  conditions.  On  the  other 
hand,  given  the  existence  of  such  a diagram,  one 
may  try  to  exploit  the  analogy  with  a further  phase 
diagram  and  ask  questions  about  the  nature  of 
transitions  between  the  different  morphologies. 


2.  Formulation  of  the  problem 

We  are  interested  in  a non-equilibrium  situation  — 
growth  of  a stable  phase  from  a metastable  one.  To 
be  specific,  we  consider  the  two-dimensional  growth 
of  a pure  substance  from  its  undercooled  melt,  where 
the  growth  is  controlled  by  the  diffusion  of  the  latent 
heat  of  freezing.  It  obeys  the  diffusion  equation  and 
appropriate  boundary  conditions  at  the  moving  (not 
known  in  advance)  interface 


u,,=Dn-(Vt/s|int-Vt/L|J  (2) 


U\int  = A-d(0)K  (3) 

The  indices  ‘L’  and  ‘S’  refer  to  the  liquid  and  solid, 
respectively.  The  specific  heat,  cp,  and  the  thermal 
diffusion  constant  D are  considered  to  be  the  same  in 
both  phases,  L is  the  latent  heat;  U = (T  — T.f)cp/L 
is  the  appropriately  rescaled  temperature  field  mea- 
sured from  the  imposed  temperature  T„  of  the 
undercooled  melt  far  away  from  the  interface;  in 
terms  of  these  parameters, 

A={TM-T„)cp/L  (4) 

is  the  dimensionless  undercooling  of  the  melt  and  TM 
is  the  melting  temperature. 

The  physics  underlying  Eqs.  (l)-(3)  is  quite 
simple.  A solidifying  front  releases  latent  heat  which 
diffuses  away  as  expressed  by  (1);  requiring  heat 
conservation  at  the  interface  gives  (2)  (n  is  the 
normal  to  the  interface).  Eq.  (3)  is  the  local  equilib- 
rium condition  at  the  interface  which  takes  into 
account  the  Gibbs-Thomson  correction;  K is  the 
two-dimensional  curvature  and  d(0)  is  the  so-called 
anisotropic  capillary  length  with  an  assumed  four- 
fold symmetry, 

d{0)  = d0(l  - e cos  46>)  (5) 

Here  d0  = yTMcp/L 2 is  a capillary  length  propor- 
tional to  the  isotropic  part  of  the  surface  energy  y:  0 
is  the  angle  between  the  normal  n to  the  interface 
and  some  fixed  crystallographic  direction,  at  which 
d(0)  is  minimal;  e is  the  strength  of  the  anisotropy. 

In  Eq.  (3)  we  neglect  the  kinetic  effects,  that  is  the 
dependence  of  the  interface  temperature  on  the 
growth  velocity  vn  which  holds  at  the  sufficiently 
small  undercoolings  and  velocities. 

Our  main  interest  here  is  concerned  with  patterns 
which  can  grow  at  constant  speed  even  at  low 
undercoolings  A<  1,  because  if  they  exist  they  will 
dominate  the  systems  behavior.  A two-phase  struc- 
ture then  must  exist  behind  the  growth  front  filling 
the  space  uniformly  on  sufficiently  large  scales.  The 
fraction  77  of  solid  inside  this  two-phase  region 
should  be  equal  to  A due  to  global  conservation, 

V = A (6) 

One  may  define  an  envelope  over  the  front  of  this 
complex  two-phase  structure,  calling  this  suitably 
averaged  envelope  the  average  front  in  contrast  to 
the  local  interface  separating  the  solid  from  the 
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liquid.  This  average  front  can  be  considered  as  the 
real  growth  front  in  the  sense  that  a two-phase 
mixture,  solid  plus  liquid,  grows  into  a one-phase 
region  originally  consisting  of  liquid  only.  These 
two-phase  structures  are  developed  from  initially 
smooth  interfaces  by  the  well-known  Mullins- 
Sekerka  instability  [15]. 

Eqs.  ( 1 )— (5)  contain  two-dimensional  parameters, 
d0  and  D,  and  two  dimensionless  parameters,  A and 
e.  It  means  that  any  characteristic  length  scale, 
and  growth  velocity,  v , of  the  possible  structures  can 
be  presented  in  the  form 

£ = d0f(A , e),v=-r  <p(A,  e ) (7) 

“o 

Our  aim  is  to  predict,  for  given  undercooling  A and 
anisotropy  e,  the  type  of  the  two-phase  structure,  and 
its  characteristic  length  scales  and  velocity,  that  is  to 
calculate  the  functions/ and  cp  in  the  relation  (7).  As 
it  turns  out  these  functions  have  scaling  forms  for 
small  A and  e,  thus  showing  power  law  dependencies 
on  A and  e. 

We  construct  the  kinetic  phase  diagram  in  the 
plane  (A,  e)  (Fig.  2),  which  represents  the  regions  of 
existence  of  different  structures  and  the  lines  of 
transitions  between  the  structures  [1,2].  We  discrimi- 


Fig. 2.  Kinetic  phase  diagram. 


nate  between  compact  structures  (C)  and  fractal 
structures  (F).  A complementary  classification  deals 
with  the  existence  of  orientational  order.  A structure 
with  pronounced  orientational  order  will  be  called 
dendritic  (D),  and  without  apparent  orientational 
order  it  will  be  called  seaweed  (S). 

It  turns  out  that  noise  which  always  exists  in  the 
system  (for  example  the  thermodynamic  noise)  ap- 
pears to  play  a crucial  role  in  the  formation  of  fractal 
structures  but  is  not  so  important  for  compact 
patterns. 


3.  Compact  dendrites  (CD) 


Dendrites  can  grow  at  arbitrary  small  undercooling 
A,  but  usually  a non-zero  value  of  the  anisotropy  e is 
required.  The  growth  pattern  evolving  from  a nu- 
cleus acquires  a star-shaped  envelope  surrounding  a 
well-defined  backbone.  The  distances  between  the 
comers  of  the  envelope  increase  with  time.  For  small 
undercooling  we  can  use  the  scaling  relation  for  the 
motion  of  the  comers  as  for  free  dendrites  [16-19] 
with  tip  radius  p,  and  velocity  v.  These  two  relations 
come  from  the  Ivantsov  formula  [20] 


(8) 


and  from  the  selection  condition  for  the  stability 
parameter  cr. 


I _ P\  ~ e~1,H 

'fir  \/d0Dlv  € 


(9) 


It  is  quite  remarkable  that  Eq.  (8)  was  obtained  in 
1947  but  it  took  about  40  years  to  derive  the  very 
non-trivial  relation  (9).  One  can  find  the  details  in 
[16-19].  From  Eqs.  (8)  and  (9)  follow  the  dependen- 
cies of  p,  and  v on  the  parameters  A and  e: 


P,~doe 


A~-,v 


1,4  A4 


(10) 


Eq.  (10)  really  describes  a needle-crystal  which, 
without  noise,  has  no  sidebranches.  The  corre- 
sponding star  structure  then  cannot  fill  the  space  with 
constant  density  and  the  amount  of  material 
solidified  in  parabolic  form  increases  with  time  only 
according  to  t:n  rather  than  t2  for  a truly  compact 
object  in  two  dimensions. 
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A small  amount  of  noise,  however,  cures  this 
problem.  The  tip  of  the  dendrite  is  still  stable  against 
small  noise  but  has  a 'convective’  instability  which 
produces  sidebranches.  Those  branches  continue  to 
grow  until  they  become  independent  primary  branch- 
es a distance  £ = Dlv  away  from  the  comers  of  the 
star.  The  global  shape  then  consists  of  an  envelope 
of  diamond  type  over  the  dendrite  tips  which  appear 
a distance  Dlv  apart  from  each  other.  The  velocity 
scales  like  Eq.  (10).  The  relative  space  filling  by 
primary  dendrites  and  sidebranches  of  course  must 
be  equal  to  A.  The  two  basic  length  scales  in  this 
pattern  accordingly  are  the  diffusion  length  Dlv  and 
the  tip  radius  pt  of  a typical  dendrite. 

For  small  undercooling  A those  two  length  scales 
are  well  separated,  pt  « Dlv.  While  the  dendritic 
structure  becomes  compact  only  at  length  scales 
larger  than  Dlv,  it  shows  fractal  behavior  at  the 
intermediate  length  scale  £,  pt<  £ <D/v,  with 
fractal  dimension  Df  = 1.5  [24].  In  this  fractal  object 
the  sidebranches  interact  due  to  the  competition  in 
the  common  diffusion  field.  Some  of  the  sidebran- 
ches die  and  some  continue  to  grow  in  the  direction 
prescribed  by  the  anisotropy.  This  competition  leads 
to  coarsening  of  the  stmcture  in  such  a way  that  the 
distance  between  the  surviving  sidebranches  is  ad- 
justed to  be  of  the  same  order  of  magnitude  as  the 
length  of  the  sidebranches  and  is  proportional  to  the 
distance  from  the  dendritic  tip.  At  the  same  time,  the 
thickness  of  the  surviving  sidebranches  is  propor- 
tional to  the  square  root  of  the  product  of  p,  and  the 
distance  from  the  tip.  On  length  scales  larger  than 
Dlv  the  dendritic  structure  appears  to  be  compact 
with  mean  density  17  = A. 


4.  Compact  seaweed  (CS) 

CD  structures  formally  exist  at  arbitrary  small 
anisotropy  e but  their  velocity  goes  to  zero  as  e 0. 
It  was  recently  discovered  that  there  is  another 
structure,  compact  seaweed  (CS),  which  is  favorable 
for  smaller  e and  larger  A.  The  velocity  of  the 
structure  remains  finite  at  e = 0. 

The  compact -seaweed  morphology  [1,2]  was 
originally  introduced  on  the  basis  of  experimental 
observations  under  the  name  dense  branching  mor- 
phology [25].  At  that  time,  however,  its  introduction 


as  a morphological  ‘phase’  distinct  from  the  well- 
known  dendritic  morphologies  was  rather  specula- 
tive. Computer  simulations  also  were  inconclusive  at 
that  time. 

The  first  indication  for  the  existence  of  such  a 
distinct  phase  came  (to  our  knowledge)  from  argu- 
ments [1,2]  based  on  a theoretical  study  of  crystal 
growth  in  a channel  [26].  This  analysis  of  channel 
growth  gave,  among  other  things,  the  following 
results.  A finger  type  pattern  symmetrically  in  the 
center  of  the  channel  could  grow  at  a constant 
growth  rate  for  dimensionless  supercoolings  A > 0.5. 
The  finger  looks  similar  to  the  Saffmann-Taylor 
finger  of  viscous  flow,  but  belongs  to  a different 
branch  of  the  mathematical  solution.  The  growth  rate 
of  the  crystal  increases  with  increasing  driving  force 
A,  as  expected.  A specifically  remarkable  result  of 
this  theory  [26]  is  that  the  driving  force  sets  a length 
scale  and  thereby  also  a velocity:  For  a given  driving 
force  0.5  < A < 1 there  exists  a characteristic  chan- 
nel width  below  which  such  a steadily  growing 
finger  is  no  longer  possible. 

However,  it  has  been  discovered  recently  that  the 
spectrum  of  solutions  for  growth  in  a channel  is 
much  richer  than  had  previously  been  assumed. 
Parity-broken  solutions  were  found  [27]  and  studied 
numerically  in  detail  [28-30].  A similar  solution 
exists  also  in  an  infinite  space  which  was  called 
‘doublon’  for  obvious  reasons  [28,29].  It  consists  of 
two  fingers  with  a liquid  channel  along  the  axis  of 
the  symmetry  between  them.  It  has  a parabolic 
envelope  with  radius  pl  and  a liquid  channel  of 
thickness  h.  The  Peclet  number,  P = vpJID  depends 
on  A according  to  the  Ivantsov  relation  (8).  The 
analytical  solution  of  the  selection  problem  for 
doublons  [31]  shows  that  this  solution  for  isotropic 
systems  (e  = 0)  exists  even  at  arbitrary  small  under- 
cooling A and  obeys  the  following  selection  con- 
ditions: 


fa 


Pt 


\/d0D/v 


Eqs.  (8)  and  (11)  give 

-7  Do 

pt~h~  d0A  , v—j-A 


(ID 


(12) 


28 


E.  Brener  et  al.  / Solid  State  Ionics  131  (2000)  23-33 


If  one  includes  finite  anisotropy  e,  doublon  solutions 
exist  only  above  the  solid  line  on  Fig.  2,  for  which 

4~e fIM  (13) 

For  e smaller  than  that  given  by  (13)  the  doublons 
obey  the  same  scaling  law  as  given  by  Eq.  (12)  [31]. 

It  should  be  noted  that  doublons  in  the  range  of 
their  existence  (13),  grow  faster  than  dendrites  for 
the  same  parameters  A and  e.  This  statement  is 
confirmed  by  numerical  calculations  [28,29], 

The  numerical  calculations  also  show  that  the 
double-fingering  structure  is  stable  against  competi- 
tion between  the  two  fingers  which  belong  to  the 
doublon.  It  means  that  the  axis  of  symmetry  and  the 
direction  of  growth  are  stable.  Of  course  these 
directions  are  arbitrary  in  isotropic  systems.  It  is  not 
completely  clear  at  the  moment  if  the  stability  of  the 
free  doublon  pair  follows  precisely  the  scaling  law 
Eq.  (13).  In  any  case  this  line  represents  a lower 
bound  on  A for  a given  e. 

We  assume  that  the  doublons  seem  to  represent  a 
key  point  in  the  growth  of  compact-seaweed  mor- 
phology (Fig.  3).  The  formation  of  a full  CS-struc- 
ture  evolving  from  a growing  nucleus  is  possible 
only  due  to  noise,  which  triggers  sidebranches,  as  in 


Fig.  3.  Seaweed  structure. 


the  CD-structure.  The  resulting  two-phase  structure 
has  an  almost  isotropic  circular  envelope  which 
moves  with  approximately  the  same  velocity  [Eq. 
(12)]  as  a free  doublon.  The  structure  is  fractal  with 
D f =1.5  in  the  intermediate  length  scale  between  p, 
and  D/v,  and  it  becomes  compact  with  mean  density 
r)  = A for  the  length  scale  larger  than  D/v.  The 
region  above  the  solid  line,  A~eU4  (Fig.  2),  corre- 
sponds to  CS-structure  where  doublons  exist  and 
grow  faster  than  dendrites.  This  line  represents  the 
discontinuous  transitions  between  CD  and  CS-struc- 
tures  with  a jump  of  velocities. 

5.  Surface  dewetting  as  a diffusional  growth 
process 

As  an  application  of  this  theory  the  effective 
equations  of  motion  for  a drying  thin  film  wetting  a 
substrate  are  derived.  These  equations  are  equivalent 
to  the  one-sided  model  of  diffusional  growth  with  an 
effective  diffusion  coefficient  which  depends  on  the 
viscosity  and  on  the  thermodynamic  properties  of  the 
thin  film. 

According  to  the  description  given  in  [33],  and 
also  by  Sharma  [34,35]  and  de  Gennes  [36],  there  is 
a possibility  for  the  almost  dry  part  of  the  solid 
substrate  to  be  in  equilibrium  with  the  wet  part 
which  is  in  fact  a thin  (but  macroscopic)  film  of  a 
liquid.  Both  parts  (dry  and  wet)  on  the  solid  substrate 
are  separated  by  an  interface,  which  can  be  described 
by  a height  variable  h( x)  with  x being  the  coordinate 
across  the  interface  from  the  dry  to  the  wet  part. 
Towards  the  dry  part,  the  height  variable  goes  to  a 
very  small  value  h , towards  the  wet  part  the  film 
thickness  goes  to  an  equilibrium  value  hx(p)  for 
given  pressure  p in  coexistence  with  the  vapor  phase. 
At  a specific  pressure  p0  the  liquid  film  can  be  also 
in  equilibrium  with  the  (almost)  dry  surface,  the 
corresponding  thickness  of  this  wet  film  then  is 
defined  as  h(l  = h,{ptA-  For  lower  vapor  pressure, 
the  equilibrium  film  thickness  h,(p)  would  decrease 
to  a value  smaller  than  h„,  but  it  would  be  metastable 
only  and  the  the  stable  dry  area  would  expand  at  the 
cost  of  the  wet  area.  This  is  the  dewetting  phenom- 
enon under  consideration.  Patterns  similar  to  Fig.  3 
have  been  recently  observed  experimentally  [37]. 

We  assume  a surface-tension  y exists  between  the 
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liquid  and  the  vapor.  The  free  energy  of  the  film  then 
can  be  written  as 

G = J {g(h(x,  y))  + ~ \W)  dr  dy  (14) 

In  equilibrium  a double-tangent  construction  to  g(h ) 
gives  two  solutions  of  h_  for  the  dry  part  and  of 
h0(p)  for  the  wet  part  of  the  surface.  This  leads  to 
the  evolution  equation  [34,35,37]  for  the  film  vari- 
able h(x,y;t): 


Note  that  a relaxational  term  proportional  to  a has 
been  added.  This  term  alone  guarantees  that  a 
homogeneous  liquid  film  will  relax  to  its  equilibrium 
value  by  evaporation  or  condensation.  /j,(p)  is  the 
chemical  potential  of  the  vapor.  For  h = hjp)  this 
term  vanishes. 

The  first  part  of  Eq.  (15),  proportional  to  the 
inverse  viscosity  of  the  liquid  film,  describes  a 
creeping  motion  of  a thin  film  flow  on  the  surface.  In 
the  (almost)  dry  area  the  contributions  of  both  terms 
to  the  total  flow  and  evaporation  of  material  can  be 
basically  neglected,  because  of  the  small  value  of 
h_,  typically  less  than  one  monolayer  of  adsorbed 
fluid.  Inside  the  wet  area  we  can  to  the  lowest  order 
linearize  h = h^[  1 + u(x,  y)],  where  u is  now  a small 
deviation  from  the  asymptotic  equilibrium  value  for 
h„(p)  in  the  liquid.  Since  Vfioo(p)  = 0 the  only 
surviving  terms  are  linear  in  u and  its  spatial 
derivatives  Vn  and  Au.  Therefore,  inside  the  wet  area, 
the  evolution  equation  for  the  variable  part  u of  the 
height  variable  h becomes 

du 

~ Deff  An  — Aeffn  (16) 

We  have  dropped  here  the  terms  ~yA 2u  since  the 
effective  diffusion  constant  Deff  = Qi\l?>r}){d1gl 
d/?2}  + ay  is  positive  and  dominates  the  long-wave- 
length behaviour  over  the  fourth-order  term.  The 
relaxation  coefficient  is  Aeff  = a{d2g/d/r2}.  Deriva- 
tives are  taken  around  the  equilibrium  value  h = h0. 
Note  that  (16)  now  is  precisely  the  equation  of 


motion  studied  in  [41],  generalizing  our  basic  model 
(Eq.  (1)). 

Of  course,  this  approximation  holds  only  inside 
the  wet  region,  not  directly  at  the  dry-wet  interface. 
This  interface  region  gives  rise  to  a profile  h(x) 
similar  to  a tanh  function.  The  development  of  the 
tanh  profile  from  the  dry  to  the  wet  part  occurs  over 
a distance  that  is  short  compared  to  the  typical 
patterns  being  observed  in  the  dewetting  process.  We 
can  therefore  replace  this  profile  by  a sharp  interface 
between  the  dry  and  the  wet  part,  but  must  add  the 
corresponding  boundary  conditions  to  the  equation  of 
motion  (16)  for  the  wet  side.  Obviously,  the  bound- 
ary conditions  consist  of  a conservation  law  Eq.  (2) 
which  guarantees  that  a displacement  of  the  dry-wet 
interface  must  locally  conserve  the  fluid.  Under 
dewetting  conditions  this  leads  to  a swelling  of  fluid 
u > 0 at  the  interface.  The  second  condition  clearly 
comes  from  the  surface  tension  y which  tends  to 
keep  the  dry-wet  interface  straight.  This  is  then  just 
the  usual  Gibbs-Thomson  condition  for  an  interface, 
as  described  in  Eqs.  (3)  and  (4)  above,  with  the 
capillary  length  being  approximately  d0  ~ yl 
({d2g(h)/dh2}£)  where  £ is  the  thickness  of  the 
dry-wet  interface,  and  the  dimensionless  driving 
force  is  A — (h0  - h„)/h0.  In  summary,  we  have  for 
this  viscous  fluid-flow  problem  of  surface  dewetting 
exactly  the  same  equations  as  for  the  diffusional 
growth  of  an  isotropic  solid. 


6.  Fractal  structures 


For  the  compact  structures  described  above  noise 
is  important  only  as  the  trigger  of  sidebranches.  It 
has  been  assumed  that  the  tips  (of  dendrites  or 
doublons)  remain  undestroyed.  However,  the 
strength  of  noise  may  be  large  enough  not  only  to 
trigger  the  sidebranches  but  also  to  destroy  the  tips. 
In  order  to  estimate  the  parameters  for  which  it 
happens  let  us  look  at  the  theory  of  sidebranch 
formation  more  carefully.  According  to  the  result  of 
Langer  [32]  the  root  mean  square  amplitude,  {£')' / > 
of  the  sidebranches  on  the  underlying  parabolic 
interface  generated  by  thermal  fluctuations  depends 
on  the  distance  from  the  tip  z according  to 
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Here  the  stability  parameter  a is  given  by  Eqs.  (9) 
and  (11)  for  dendrites  and  doublons,  respectively;  r 
is  the  relative  noise  strength  if  «:  1) 


r = (T/T0) 


2 Dell 
__ 

vp  t 


T = 
» 1 o 


L-dj 

Vp 


(18) 


where  kB  is  the  Boltzmann  constant.  The  tip  be- 
comes destroyed  if  the  amplitude  of  the  sidebranches 
is  of  the  order  of  pt  at  the  distance  z ~ p,  down  the 
shaft.  Thus  we  obtain  from  Eq.  (14)  the  following 
condition 


(19) 

The  tips  of  the  structures  will  be  destroyed  if  the 
stability  parameter  a becomes  smaller  than  the 
critical  value  cr*  given  by  Eq.  (16).  Using  the  value 
cr~e7/4  [Eq.  (9)],  one  obtains  from  Eq.  (16)  a line 
of  smooth  transition  from  CD  to  FD-structures  in 

Fig.  2: 


lins-Sekerka  instability.  At  that  time,  however,  the 
existence  of  doublons  was  not  known. 

A new  approach,  therefore,  is  required  for  the 
description  of  the  fractal  patterns  with  the  destroyed 
tips.  Such  destroyed  fractal  structures  have  been 
already  investigated  in  the  framework  of  Saffmann- 
Taylor  viscous  fingering  and  diffusion-limited  aggre- 
gation [21-23],  The  important  result  of  these  in- 
vestigations is  that  there  exists  an  effective  envelope 
obtained  by  averaging  over  the  structures,  which  has 
precisely  the  same  shape  as  an  ideal  stable  solution 
— the  shape  of  the  Saffmann-Taylor  finger  in 
isotropic  systems  and  parabolic  shape  in  anisotropic 
systems.  The  density  inside  this  effective  envelope  is 
f]  < 1 . The  envelope  has  a characteristic  tip  radius 
pr  Because  the  underlying  structure  is  fractal  with 
fractal  dimension  £>,  = 1.71  in  the  intermediate 
length  scale  between  small-length  cutoff  ar  and  p„ 
the  density  f)  inside  the  envelope  can  be  obtained 
from  the  definition  of  the  fractal  dimension  (apart 
from  a constant  prefactor) 


e*  ~ |ln  r\~*n  (20) 

The  analogous  line  which  separates  CS  and  FS- 
structures  in  Fig.  2 can  be  obtained  using  Eqs.  (8), 
(11)  and  (16): 


t\ 

d r rfj(r) 

"/■ 

which  gives  more  explicitly 


A*~\ln  r\~215  (21) 

Eq.  (16)  has  the  following  physical  meaning.  Let  us 
rewrite  this  relation,  using  the  definition  of  cr  = d0D/ 
(vp f)  [Eq.  (9)],  which  gives  the  following  condition 
for  a stable  tip  radius: 

Pt-PMslln/1  (22) 

where  pMS  ~\/df)Dlv  is  the  Mullins-Sekerka  length 

describing  the  instability  of  a planar  interface.  One 
can  think  of  the  right-hand  side  of  (19)  as  the 
characteristic  length  scale,  ar,  of  the  instability  due 
to  noise 

aT  = Pms  lln  ~\jd{)D/v  |ln  r\  (23) 

The  tip  is  stable  if  pt<ar  and  it  becomes  destroyed 
if  pt  > ar.  The  same  small-size  cutoff  ar  depending 
on  noise  r was  introduced  in  our  previous  paper 
[1,2]  based  on  the  consideration  of  the  self-similar 
development  of  perturbations  induced  by  the  Mul- 


(24) 


Following  these  results,  we  will  now  try  to  define  an 
averaged  or  coarse-grained  structure  over  such  a 
noisy  fractal  pattern  and  to  formulate  an  equation  of 
motion  for  this  coarse-grained  structure  using  scaling 
arguments.  More  explicitly  we  try  to  estimate  the 
characteristic  length  scale  of  the  structure  and  its 
growth  velocity  by  considering  the  steady-state 
motion  of  an  effective  parabolic  envelope  which 
replaces  the  destroyed  dendrite  or  doublon.  The 
density  of  the  solid  phase  inside  the  envelope  with 
tip  radius  px  is  given  by  Eq.  (21),  where  the  small 
size  cutoff  ar  is  defined  by  Eq.  (20) 


77(P,)~ 


/ y/d0D/v\]n  r\ 
V Px 


2-0, 


(25) 


The  temperature  inside  the  envelope  is  assumed  to  be 
close  to  the  melting  temperature.  Because  the  density 
fj  inside  the  envelope  is  smaller  than  1 , we  have  to 
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replace  the  latent  heat  L by  pL.  It  changes  A in  Eq. 
(4)  to  {Alp)  and  modifies  the  Ivantsov  relation  to 


P 


(AY 

2 D~\r,) 


(26) 


The  crucial  point  of  the  analysis  is  a modification  of 
the  selection  conditions  (9)  and  (11).  The  ex- 
perimental and  numerical  results  [21-23]  (the  exist- 
ence of  a selected  envelope)  support  the  idea  that 
those  selection  conditions  do  exist.  Unfortunately, 
we  do  not  know  any  results  which  allow  us  to  write 
down  these  modifications  explicitly.  But,  using  sca- 
ling arguments,  we  can  write  the  selection  conditions 
in  the  following  scaling  form  with  scaling  exponent 
(3  which  for  the  moment  is  undetermined  but  will  be 
specified  later.  For  dendrites  we  can  write 


1 _ Pt -7/8  ~-/3 

■Jo-  ~ \Jd0Dlv  ~ € V 

and  for  doublons 

J-  = K „p-si 4 

'fc  \Jd0D/v 


(27) 


(28) 


These  relations  (24)  and  (25)  transform  into  Eqs.  (9) 
and  (11),  respectively,  for  p~l. 

The  selection  relations  (24)  and  (25)  may  be 
interpreted  as  the  conditions  of  selection  due  to  an 
effective  surface  tension.  We  have  chosen  the  same  (3 
in  both  Eqs.  (24)  and  (25)  because  in  some  sense  the 
factor  p ^ can  be  seen  as  a renormalization  factor 
for  the  capillary  length  da  —>  d0  / p >fl.  To  make  an 
estimate  of  possible  values  of  f3  it  is  natural  to 
assume  that  a coarse-grained  surface  energy  should 
decrease  with  decreasing  p,  giving  /3  < 1/2  as  a 
reasonable  restriction.  We  have  shown  [1,2]  that  the 
scaling  exponent  (3  can  be  expressed  in  terms  of  bulk 
Dt  and  surface  Ds  fractal  dimensions: 


A-l 

2 — D{ 


(29) 


All  the  characteristics  of  fractal  structures  depend  on 
the  noise  strength  F We  can  estimate  r in  the  fractal 
region  using  Eq.  (15)  and  replacing  px  by  the  noise- 
induced  length  scale  ar  from  Eq.  (20).  It  gives 


r|in  r |2~ 


1/2 


(30) 


The  noise  strength  r •«  1 because  the  capillary 
length  da  is  much  smaller  the  diffusion  length  Dlv. 


7.  Conclusion 

We  have  discussed  the  structure  formation  in 
diffusion-controlled  growth.  The  given  description 
refers  to  solidification  of  a pure  undercooled  melt  but 
it  also  can  be  applied  to  growth  of  a pure  solid  from 
solution  or  isothermal  solidification  of  a binary  melt. 
More  generally  one  may  speak  of  systems  with  a 
conserved  quantity  growing  by  diffusion.  The  main 
control  parameters  of  the  process  are  dimensional 
undercooling  A and  the  strength  of  the  surface 
tension  anisotropy  e.  It  turns  out  that  the  noise  is  also 
very  important  for  the  structure  formation  and  we 
characterize  it  by  the  dimensionless  quantity  r. 

The  resulting  morphology  diagram  (Fig.  2)  with 
axes  4 vs.  f classifies  different  kinds  of  structures 
and  transitions  between  them. 

The  dendritic  structure  has  pronounced  orienta- 
tional order  and  it  is  favorable  for  small  A and 
relatively  large  e.  The  seaweed  structure  does  not 
require  anisotropy  and  is  favored  for  larger  A and 
smaller  e.  The  transition  between  these  two  struc- 
tures takes  place  around  the  solid  line  on  Fig.  2 [Eq. 
(13)]  which  is  continued  by  the  dotted  line  into  the 
fractal  region.  This  transition  is  discontinuous  with  a 
jump  of  velocities  since  the  doublons  move  faster 
than  the  dendrites  as  soon  as  they  exist.  The  main 
element  of  the  dendritic  structure  is  a dendrite  with  a 
parabolic  tip,  and  the  main  element  of  the  seaweed 
structure  is  a doublon  (Fig.  3).  For  compact  dendritic 
and  compact  seaweed  structures  the  tips  of  dendrites 
and  doublons  are  stable  against  the  noise  which  is 
relatively  small  in  these  regions.  The  noise  triggers 
sidebranches  which  fill  the  space  and  make  the 
structures  compact  so  that  the  mean  density  of  the 
solid  phase  is  77  = A on  the  length  scale  larger  than 
Dlv.  In  the  intermediate  region  of  lengths  between 
the  tip  radius  pt  and  diffusion  length  Dlv  the 
structure  can  be  described  as  a fractal  but  with  a 
trivial  fractal  dimension  D{  = 3/2  which  comes  from 
the  parabolic  shape  of  the  dendrite. 

The  region  of  fractal  dendritic  and  fractal  seaweed 
structures  near  the  origin  of  the  morphology  diagram 
is  characterized  by  noise  being  sufficiently  large  to 
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destroy  even  the  tips  of  dendrites  and  doublons.  It 
means  that  the  noise-induced  length  scale  ar  [Eq. 
(20)]  is  smaller  than  pr  In  the  range  between  a,  and 
pt  the  structures  are  fractal  with  a non-trivial  fractal 
dimension  Df  (, Dr  =1.71).  This  is  the  reason  why  we 
called  these  structures  ‘fractal’.  Furthermore,  these 
patterns  are  also  fractal  in  the  range  between  pt  and 
D/v  but  again  with  the  trivial  fractal  dimension 
Df  = 3/2.  Finally  they  become  compact  on  length 
scales  larger  than  the  diffusion  length  Dlv  just  as 
compact  dendritic  and  compact  seaweed  structures. 
Note  that  if  one  performed  a measurement  of  the 
fractal  dimension  on  length  scales  around  the  cross- 
over length  p,  one  would  observe  an  interpolation 
between  our  two  different  values  of  D(,  the  precise 
result  depending  on  the  interval  chosen  for  the 
measurement.  Since  both  the  dendritic  and  the 
seaweed  patterns  maintain  their  basic  identities  in- 
side the  noisy  region  the  transitions  from  the  com- 
pact to  the  fractal  regions  represent  rather  smooth 
changes  in  length  scales. 

We  have  described  the  structures  and  growth 
velocities  of  the  destroyed  fractal  dendrites  and 
doublons  by  introducing  renormalized  quantities  for 
capillary  length  and  density.  We  have  quantitatively 
introduced  an  effective  parabolic  envelope  following 
the  results  of  [22-24].  The  most  non-trivial  part  of 
our  analysis  is  a modification  of  the  selection 
conditions  [Eqs.  (24)  and  (25)].  At  this  point  we 
have  used  scaling  arguments  which  leave  us  only 
with  one  undetermined  scaling  exponent  fi.  This 
exponent  subsequently  is  determined  by  the  fractal 
dimensions  Z)s  for  the  surface  and  D,  for  the  bulk  of 
the  growing  pattern.  The  closed  set  of  equations  for 
the  growth  rate  and  the  tip  radius  are  Eqs.  (22)-(24) 
for  fractal  dendritic  growth,  and  Eqs.  (22),  (23)  and 
(25)  for  fractal  doublon  growth.  The  two-dimension- 
al theory  together  with  a small  modification  [41]  in 
addition  explains  the  patterns  observed  on  a solid 
surface  which  was  covered  by  a thin  wetting  layer  of 
fluid.  When  a dry  spot  nucleates  in  this  layer 
dewetting  occurs  through  the  propagation  of  the 
essentially  one-dimensional  separation  line  between 
the  wet  and  the  dry  parts  of  the  surface.  Doublon 
structures  then  can  be  clearly  observed. 

The  scaling  arguments  given  here  for  two-dimen- 
sional growth  patterns  formally  can  be  extended  in  a 
straightforward  fashion  to  three  dimensions.  For 


dendritic  structures  this  seems  to  be  perfectly  permis- 
sible since  the  basic  growth  laws  are  rather  similar  in 
two  and  three  dimensions  [38,39]: 

p{~d0  e~"A  |ln  A\IA,  v ~~r  e7/4  (A/\n  zl)2  (31) 
“a 

(compare  to  Eq.  (10)  for  2-D).  There  is,  however,  the 
crucial  difference  between  the  3-D  and  the  2-D  case. 
In  the  later,  small  anisotropy  implies  that  the  shape 
of  the  selected  needle  crystal  is  close  to  the  Ivantsov 
parabola  everywhere;  in  the  former,  strong  devia- 
tions from  the  Ivantsov  paraboloid  appear  for  any 
anisotropy.  This  shape,  in  units  of  the  tip  radius  of 
curvature,  depends  mostly  on  the  crystalline  symme- 
try and  it  is  almost  independent  of  the  material  and 
growth  parameters.  The  shape  of  the  3-D  dendrite, 
which  has  been  described  analytically  [39],  together 
with  sidebranching  activity  [24]  is  presented  in  Fig. 
4.  For  the  seaweed  patterns  much  less  is  known  since 
our  preliminary  results  are  mostly  numerical  ones 
[40]  (Fig.  5).  The  crucial  point  here  is  that  this 
self-organized  triplet  structure  is  not  imposed  by  the 
symmetry  of  the  calculation  box  and  it  consists  of 
three  cooperating  symmetry-broken  fingertips.  A 
hexagonal  or  triplet  structure  should  be  expected  to 
occur  under  free  growth  conditions  from  the  basic 
symmetry  considerations,  since  these  growth  prob- 
lems do  not  have  reflection  symmetry  about  some 


Fig.  4.  3-D  dendrite. 
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Fig.  5.  3-D  triplon. 


average  interface  position.  The  resulting  multiplets 
consisting  of  symmetry-broken  cooperating  fingers 
of  the  growing  phase  finally  seem  to  be  the  basic 
building  blocks  for  the  compact  seaweed  morpholo- 
gy, in  particular  the  triplet  structure  appears  to  be 
characteristic  for  free  growth.  We  expect  these 
triplons  even  to  appear  in  systems  with  convection 
driven  by  surface  tension,  provided  the  growth  rate  is 
faster  than  the  convection  rate  defined  as  the  ratio  of 
surface  tension  divided  by  viscosity.  At  present, 
however,  there  is  still  a need  for  detailed  experi- 
ments in  this  region  of  low  crystalline  anisotropy  but 
rather  high  driving  force. 
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